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We study flow phenomena in relativistic heavy-ion collisions, both 



in transverse and radial direction, in comparison to experimental data. 
The collective dynamics of the nucleus-nucleus collision is described 
within a transport model of the coupled channel BUU type (CBUU). 
This recently developed version includes all nucleonic resonances up to 
1.95 GeV in mass and mean-field potentials both of the Skyrme and 
momentum dependent MDYI type. We find that heavy resonances 
play an important role in the description of transverse flow above 
1 AGeV incident energy. For radial flow we analyse reaction times and 
equilibration and extract the parameters T and (3 for temperature and 
collective flow velocity within different prescriptions. Furthermore, we 
apply a coalescence model for fragment production and check the mass 
dependence of the flow signals. 
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1 Introduction 

Relativistic heavy-ion collisions (HIC) provide a unique tool to study nuclear 
matter at high densities and temperatures. However, throughout a collision 
the system is partly far away from equilibrium and both particle production 
and collective motion depend on various quantities, such as the stiffness of 
the equation-of-state (EoS), the momentum dependence of the interaction 
or mean-field potentials (MDI), in- medium modifications of the NN cross 
section a^N and even the initial momentum distribution of the nucleons P|- 



TTf . Since these dependencies are in general not very strong or unique and 
also model parameters influence the results, it is necessary not to focus on 
a single observable alone, but to investigate the dynamical evolution of the 
HIC within a single model that is able to describe all relevant quantities. 

In this work we investigate the collective behavior of nuclear matter in 
a heavy- ion collision in the energy range from 150 AMeV to 2 AGeV for 
various systems using the CBUU model which is briefly presented in Section 
D- This model was shown to describe pion production [12|], photoproduction 
and -absorption |13], [TJJ] as well as pion- induced reactions [HJ. Recently we 



investigated isospin equilibration and its dependence on the EoS as well as 



on medium modifications of the NN cross sections |16| in comparison to 
experiments currently performed at GSI [|17[ . 

The paper is organized as follows: In Section || we describe the CBUU 
model used for the calculations where we especially focus on mean field po- 
tentials. Section |3] presents the results and discussions on transverse flow, 
while Section |] is assigned to the results on radial flow. The mass depen- 
dence of the flow quantities is investigated in Section |, whereas a summary 
follows in Section || 



2 The CBUU-Model 



2.1 Basic equations and collision term 

For our study we use the CBUU transport model which has already been in- 
troduced in an earlier publication [0 . Whereas in [ 12] we have concentrated 
on particle production and therefore on the collision term, we will focus here 
on the mean-field potential, which governs the bulk properties of the nuclear 
matter throughout a heavy-ion collision. 



The basic equation of a BUU transport model reads |T8|, [L9[ [21], |22] : 



dfi 
dt 



Pi 



+ f + V P U x {r, p[) | V r fi - V r fMf, pl)V„/i 
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= E J d3 P2 J d3 P3 J d ^ ^ {Pi +P2-P3- Pi) 

X (^34 ,^ h JA J\ h - Vl-2 fl f 2 h H) , (1) 

where fi stands for fi(r,Pi,t) denoting the single-particle phase-space distri- 
bution function for the nucleons; fi = (1 — f\) are the Pauli-blocking factors 
for fermions. Since we go beyond the pion production threshold similar 
equations have to be solved for the nucleon resonances and for the mesons 
included. These transport equations are coupled via the collision term. We, 
therefore, refer to the model as the coupled channel BUU model (CBUU). 
Note, that in Eq. (JJ) all phase-space distributions / appear with equal space- 
time argument, since the collision term is assumed to be local in space and 
time. 

As already described in |12] we include 14 nucleon resonances up to masses 
of 1.95 GeV/c 2 , i.e. the A(1232), iV(1440), JV(1520), iV(1535), A(1600), 
A(1620), JV(1650), A(1675), JV(1680), A(1700), JV(1720), A(1905), A(1910) 
and A(1950), where the resonance properties are adopted from the PDG |23| . 
The mesons incorporated are ir, rj, p and a, where the cx-meson is introduced 
to describe correlated pion-pairs with total spin J = 0. For the baryons as 
well as for the mesons all isospin degrees of freedom are treated explicitly. 

The r.h.s. of Eq. ([!]), i.e. the collision integral, describes the changes of 
fi(r,pl, t) due to two-body collisions among the hadrons (h): hi+h^ <-> h^+h^ 
and two-body decays of baryonic and mesonic resonances (R, r) to hadrons 
and mesons (m): Rh *-> h + m, r m <-> m!+m 2 . Three body final states are 
treated as two subsequent 2-body processes. The in-medium collision rate is 
represented by V12 ■ da/dQ where da/dfl is the free differential cross section 
and v±2 is the relative velocity between the colliding hadrons hi and h% in 
their center-of-mass system. Taking elastic collisions only, t> 12 • dayi^^j d£l 
equals f 34 • da^^^/dVt] in case of inelastic collisions one has either to use 
detailed balance, or, as we do, divide the cross section up into (matrix element 
x phase space factors) and use these matrix elements for the determination 
of the backward reaction ||24|| . For the most important nucleonic resonance, 
the A(1232), we use a parametrization following the result of an OBE model 
calculation by Dimitriev and Sushkov p5 for mass- and angle-differential 
cross sections. 

In the collision integrals describing two-body decays of resonances the 
product (relative velocity x cross-section x f 2 ) has to be replaced by the 
corresponding decay rate and the proper fermion blocking factors in the final 
channel have to be introduced. The factor g in Eq. ([]]) stands for the spin 
degeneracy of the particles participating in the collision whereas Yj2 3 4 stands 
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for the sum over the isospin degrees of freedom of particles 2, 3 and 4. 

We include the following elastic and inelastic baryon-baryon, meson- 
baryon and meson-meson collisions: 



NN 
NN 
NR 
NN 
R 
R 

7V(1535) 
NN 

P 
a 



NN 
NR 
NR' 

A(1232)A(1232) 

Nn 

Nun 

A(i232)vr, iV(1440)7r, Np, Na 

Nij 

NNir 

Tin (p-wave) 
7T7T (s-wave) , 



(2) 



(3) 



where the 27r-decay of a hadronic resonance, R, is introduced as subsequent 
l7r-decays (|3]) and the NN < — ► AA process (0) is parametrized according 



to Huber and Aichelin . 



2.2 Mean-field potentials 

The l.h.s. of Eq. ([!]) represents the Vlasov-equation for hadrons moving 
in a momentum-dependent field U(f,p), where r and p stand for the spa- 
tial and momentum coordinates of the hadrons, respectively. From Dirac- 



phenomenological optical-model calculations [2y|] it is known that elas- 
tic nucleon-nucleus scattering data can only be described when using proper 
scalar and vector potentials. Since this approach has proven to be numerically 
very difficult [^] we use as starting point the non-relativistic momentum- 
dependent mean field potential proposed by Welke et al. @, [J, i.e. 

u(f, p) = a p -^ + b( e&y + 2^ Uy f{r ; pl \ 2 . ( 4) 

Po \ Po J Po J 1 + 



p-p 

A 



This ansatz in principle enables to guarantee energy conservation since it 
can be derived from a potential energy density functional. However, as an 
extension of the momentum-independent Skyrme type potentials for nuclear 



matter [20. 29] the parametrization (Eh has no manifest Lorentz properties, 



whereas the latter are required for a transport model at relativistic energies. 
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To achieve this goal we evaluate the non-relativistic mean-field potential U 
for a particle in the local rest frame (LRF) of the surrounding nuclear matter 
which is defined by the frame of reference with vanishing local vector baryon 
current, j(r,t) = 0. For this the reaction volume is divided into a grid with 
1 fm 3 cell volume and the baryon current is evaluated in each cell. In the 
LRF we then transform the non-relativistic potential U(r,p) to a scalar one, 
Us, by the identification 



\Jp 2 + m 2 + U(r, p) — \Jp 2 + (m + Us(r, p)) 2 , (5) 

thus defining a local momentum-dependent effective mass m* = m + Us for 
the particle. This effective mass is used throughout our calculations for the 
baryons. The mesons are propagated as free particles; their effective mass is 
equal to their restmass, i.e. U meson = 0. 

Due to the relativistic dispersion relation (|]), the potential Us(r, p) has 
now definite Lorentz-properties. This enables us to guarantee energy conser- 
vation in each two-body collision ( N\ + N2 — > N% + N4) as 



p{ + m*i + \Jp2 + m\ 2 = y/pg + 7713 + \jpl + ml 2 (6) 

as well as in resonance decays. 

The parameters of the Potential U we fit to match the requirements 



E 
A 



f) U 

= -16 MeV, — 
po dp 



= 0, K = 210 I 260 I 380 MeV, 
U(E = 300 MeV) = and U(p = 00) = +32 MeV . (7) 



These constraints we derive from the results of a microscopic calculation from 



Wiringa et al. pOl |3T[ with Hamiltonians that describe iViV scattering data, 



few body binding energies and nuclear matter saturation properties. In f31 
the results of the calculation based on the Urbana vl4 model plus additional 
three-body interaction (UV14+TNI) is fitted by the approximation 

U = a(p)+ m 



1 + 00 

for different densities, which is close to the functional form of Eq. (J|). Thus 
we fit Eq. (H) with (|j) for normal nuclear matter density which we take to be 
Po ~ 0.168 fm -3 . The best agreement to (§) over the whole density regime 
from 0.1 to 0.5 fm -3 is achieved using a compressibility of K ~ 230 MeV, 
however, we fit the three different compressibilities denoted in Eq. (^), K — 
210, 260 and 380 MeV, to allow for simulations testing different EoS. 
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2.3 Numerical realization 



The CBUU-equation ([!]) is solved by means of the test-particle method, where 
the phase-space distribution function f(r,p,t) is represented by a sum over 
5-functions: 



NxA 



f(r, p, t) 



N 



(5(r-n(f))x^(p-ft(t)), 



(9) 



i=l 



Here N denotes the number of testparticles per nucleon while A is the total 
number of hadrons in the two colliding nuclei. Inserting the ansatz (|9]) into 
the CBUU-equation leads to the equations of motion for the testparticles: 



dfjjt) 

dt 
dpijt) 

dt 



I + — V p U s (fi, £•(*)) 



m 



— V r U s (n, Pi{t)). 



(10) 



Thus the solution of the CBUU-equation within the testparticle method re- 
duces to a set of equations of motion for classical point particles (fL0~l). For 
the actual numerical simulation we discretize the time t into steps of typ- 
ical 0.5 fm/c and integrate the equations of motion employing a predictor- 



corrector method 32 



For the evaluation of the mean-field potential we expand the local Thomas- 
Fermi approach used for the initial momentum distribution for the testpar- 
ticles. Then the integral in Eq. (01) can be solved analytically giving B] 



U(f,p) 



7rA 



Pp + A 2 - p 2 ^ (p + p F f + A 2 M 2p F 



2pA 



arctan 



(p-Pf) +A 2 A 



P + Pf 
A 



arctan 



P~Pf 
A 



(11) 



which provides a very fast evaluation of the nucleon potential. Two-body 
collisions and resonance decays are treated as described in Ref. (12[ . 

This model has been proven to adequately describe pion spectra |]12|, pion 
multiplicities |24| and Coulomb effects on pion spectra |33 . 



2.4 Energy conservation 

One of the most important issues to check for any transport model is the 
quality of energy conservation since this influences the results on both trans- 
|29H and radial flow as well as on particle production. We evaluate the 



verse 



total energy minus the nucleon restmass 
E = + Jd 3 rW(p,r) 
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i 

( 



+ d 3 r 



N 



A p(rl 2 , B P(rY T+1) , C 



\ 



Po 



Pi 



Po 



73 73 / f(r,P)f(r,p') 
a pa p 



p-p 



NxA 

E 

i=l 



mf + pf 



+ E 



/ 2 



B 1 - 



+ 



g 



'12) 



where the kinetic energy sum runs over all N x A testparticles and the spacial 
integral / dxdydz W is taken as a sum over all grid cells J2 x ,y,z an d J2j,k runs 
over all particles in the current cell (x,y,z). The momenta Pi,p~j are taken 
relative to the LRF. 

The upper part of Fig. [1| shows this total energy per nucleon in a 1 AGeV 
central Au + Au collision in the cascade mode (i.e. neglecting the mean field), 
the middle part the results for propagating the nucleons in the mean field 
without collisions and the lower part for the 'full' calculation, i.e. propagation 
of the particles in the mean-field plus collisions. 

Whereas the collisional part (upper picture) - which we require to guar- 
antee energy conservation in the individual collisions on the 10 _5 -level - 
gives rise to a steady total energy loss about 2 %, the propagation in the 
momentum- dependent mean field (middle picture) in total conserves energy 
but shows some deviation in the high density phase between t=7.5 fm/c to 
t=25 fm/c. These deviations from the initial total energy per nucleon are 
due to a gaussian smearing algorithm for the density distributions used to 
obtain smooth potentials in coordinate and momentum space. This avoids 
steep gradients for numerical reasons which would lead to an artificial and 
unphysical acceleration of the particles. The smearing, of course, smoothes 
also steep physical density gradients, e.g. in the region of the fireball surface. 
However, taking the 'bare' density distribution for the evaluation of the po- 
tentials leads to numerically instable nuclei with the consequence that the 
results depend on initial conditions as extensively discussed in Ref. ||; with 
smearing they are stable. 

The full calculation (lower picture) shows mainly the behaviour of the 
cascade with a total energy loss of about 1.8 %, which is an excellent value 
for the following investigations. 
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3 Transverse flow 

Within the CBUU model described in the previous Section we now calculate 
the transverse flow for various systems and beam energies and analyse the 
dependence on different quantities. 

In general, in a HIC three stages can be identified: i) the initial phase, 
which is characterized by high relative momenta of the incoming nucleons, ii) 
the high compression phase, where due to the high density a lot of collisions 
happen thus driving the system towards thermal and chemical equilibrium, 
and iii) the expansion phase, where still some particle freeze-out happens, 
but the collision rate drops to zero. The phases i) and ii) are most important 
for collective flow as will be shown in the following. 



3.1 Origin of flow 

In the participant - spectator picture transverse flow has its origin in the 
deflection of the spectators at the hot and dense reaction zone, the fireball 
(see, e.g. pl[ ). However, in a transport model the flow pattern of the 
nucleons in coordinate space (Fig. |2]) for a Au + Au collision at 1 AGeV 
and impact parameter b=6 fm shows that the direction of motion of the 
spectators in the final state of the reaction, clearly characterized by a density 
close to po, is nearly unchanged. At the back of these spectators participating 
matter escaping from the fireball streams in outward direction alongside the 
spectator surface, attracted by the nuclear potential field. Since the flow F is 
defined as the slope of the transverse momentum distribution at midrapidity 



d(p x ) 



dy 



(13) 



y=yo 

it is this participating matter which gives rise to the transverse flow signal. 

This can also be seen by comparing the flow values derived from the 
transport calculation when using either all nucleons or cutting out the spec- 
tator nucleons, which can easily be identified by their collision number. Now 
defining the flow F as derivative between —0.5 < y/yo < 0.5 we obtain for a 
Au + Au collisions at b=6 fm the results displayed in Table [l]. Clearly only at 
very low incident energies (~ 150 AMeV) - when the spectators move with 
relatively low velocities - the results including or excluding spectator matter 
differ when calculating the flow. 
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E, ■ 

J -'kin 


F -it [MpVI 

r with L ivlt! V J 


F i [MpVI 


1 ^o ampV 


VJVJ 


8fi 


250 AMeV 


112 


115 


400 AMeV 


138 


139 


1 AGeV 


197 


197 


2 AGeV 


222 


221 



Table 1: Flow values (dp x /dy) for Au + collisions at b=6 fm for dif- 
ferent incident energies as calculated within the CBUU model. Center row: 
Including all particles, right row: Excluding particles with collision number 
zero. 



Based on the observation, that sidewards flow is created by the partici- 
pating matter, it becomes transparent why flow does not clearly distinguish 
between an EoS with and without momentum-dependent forces. Since the 
fireball contains the stopped matter, the relative momenta in the fireball 
besides the unordered thermal motion, are small. Only when applying addi- 
tional cuts, e.g. on high transverse momenta |], [?|, i.e. by selecting particles 
escaping early from the fireball, or selecting mainly participant or spectator 
particles by appropriate cm -Cuts |35| , a difference between the momentum- 
dependent and momentum-independent EoS can be established. 



3.2 Conservation of angular momentum 

Since flow is generated by the participating matter escaping from the fire- 
ball, where the particles move random and undergo numerous collisions, the 
inclusion of an explicit angular momentum conservation mechanism for the 
two-body collisions in a transport model is of minor importance. The conser- 
vation of angular momentum or at least the conservation of the reaction plane 
in the individual particle-particle collisions is usually neglected in transport 



models. In p9| , [36|j the effects of an explicit angular momentum conservation 
mechanism have been studied, where in [[36| a substantial effect of a reaction 
plane conservation and a systematic choice of repulsive or attractive scatter- 
ing trajectories in the individual nucleon-nucleon collisions on transverse flow 
is reported. We, therefore, investigate the influence of modifications of the 
CBUU collision term on angular momentum conservation and flow. However, 
assuming both energy and angular momentum conservation in the individual 
particle-particle collisions gives rise to some conceptional problems: These 
two quantities determine the collision time and angular distribution in a 
unique way. The latter is usually simulated in line with the differential cross 
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section measurements for free NN collisions and thus should not be changed 
in favor of an arbitrary collision geometry. The first, i.e. the time restriction, 
is in contradiction to the concept of finite time steps and requires the reloca- 



tion of the colliding particles either in coordinate or momentum space [29 



This relocation again disturbs the evolution of the phase-space density and 
gives rise to an artificial motion of the particles. To avoid these problems we, 
therefore, only checked the importance of the conservation of the reaction 
plane, that is keeping the direction of Ly constant in the individual colli- 
sions between particles i and j, where similar to [36| attractive or repulsive 
particle trajectories could be chosen. 

As starting point in Fig. |^ the total angular momentum of a semicentral 
(b=4 fm) Ni + Ni system at different incident energies is shown, where we 
have calculated 

1 N 111 

w=-tf£i*ii, u=i}%+ii+i 2 r y (14) 

iV 3=1 

with ^ 

hx = -7 5Z( r zPy — r yPz) j Ijy = ■ ■ ■ , Ijz — ■ ■ ■ j (15) 

iej 

where N denotes the number of testparticles per nucleon and A is the mass 
number of the system. The sum over i in Eq. fll5] ) runs over all particles 
(baryons and mesons) of the ensemble j. At 150 AMeV the angular momen- 
tum in the cascade mode is conserved better than 2.6 %; the mean field part 
again gives rise to some deviation from the initial value throughout the high 
density phase but in total conserves angular momentum better than 0.2 %. 
The full calculation conserves \L\ by 4.7 %. At 2 AGeV the situation looks 
similar, but the deviations are smaller. The cascade mode deviates by about 
2.2 %, the mean field mode about 0.2 % and the full BUU mode about 3.2 %. 
These rather small deviations from total angular momentum conservation al- 
ready indicate a small influence of a different treatment of the collisions on 
the evolution of the system. 

Now, including the conservation of the reaction plane in the individual 
collisions, the result is shown in Fig. f|, upper part. The quality of total angu- 
lar momentum conservation in the cascade mode is the same, 2.6 % compared 
to 2.7 %; a similar behaviour is obtained for the full BUU mode. Thus, for 
the total angular momentum the treatment of the individual collisions seems 
to play a minor role. 

In order to check a possibly more sensitive quantity, we have additionally 
calculated the azimuthal distribution of particles in the forward hemisphere, 
both with and without a conservation of the reaction plane, in the cascade 



3 TRANSVERSE FLOW 



11 



mode for semicentral collisions of systems of different mass. The results 
are shown in Fig. [|. Whereas for a pp collision the effects are expectedly 
dramatic, already for the dd system the signal is much weaker. This is due 
to the fact that the collision plane of the individual collisions is not identical 
to the reaction plane of the system and more or less arbitrary. Consequently, 
with increasing system mass and collision number the effects of a reaction 
plane conservation in the individual collisions rapidly vanishes. Table |2| gives 
the ratio R = N /N go obtained from a coscp + cos 2 (p fit to the data of Fig. ||. 



System 


b [fm] 


-Rwith 


Rw/o 


av. coll. number 


pp 


0.2 


OO 


1.02 ±0.02 


1 


dd 


1 


1.35 ±0.02 


1.00 ±0.01 


1,24 


4 He 4 He 


1 


1.22 ±0.02 


0.97 ±0.02 


1,41 


12 Q 12q 


1.5 


1.09 ±0.02 


0.96 ±0.02 


1,48 


16Q 16Q 


1.7 


1.05 ±0.01 


1.00 ±0.01 


2,35 


40 Ca 40 Ca 


2.3 


1.03 ±0.01 


0.98 ±0.01 


2,95 



Table 2: Anisotropy ratio R = Nq/N 90 at 150 AMeV for particles from the 
forward hemisphere. 



When increasing the incident energy above the one- and two-pion thresh- 
old, numerous inelastic collisions happen, where the angular momentum of 
the baryons cannot be conserved, anyway. More important than the elastic 
scattering is the treatment of the meson-nucleon-reactions here. Normally, 
when e.g. a pion matches the requirement 



' a 



tot 
ttJY 



^7riV — \ (16) 



7T 



and also fulfils the Kodama-criterion [37], it is absorbed by the nucleon and 



the created resonance is located at the position of the absorbing nucleon. Al- 
ternatively we apply a treatment, where the absorbing nucleon in a ttN — > R 
(resonance) reaction is relocated to the center of momentum of the ir and 
the nucleon. This procedure also does not conserve the angular momentum 
of the pion plus the nucleon relative to the lab system, but one can expect 
that the deviation is on average smaller than in the standard treatment. In 
Fig. |], lower part, the difference between the 'normal' BUU and this alterna- 
tive treatment is shown. The total angular momentum conservation increases 
from 2.2 % to 1.2 %. However, as in the discussion of the previous paragraph, 
this relocation gives rise to unphysical sudden motion of the baryons which 
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disturbs the time evolution of the system and which we thus do normally not 
include in our model. 

In summary, concerning sidewards flow, we find no systematic changes 
in the results of (-j^), at least in the range —0.5 < y/yo < 0.5, within the 
statistical accuracy, both for the Ni + Ni and for the Au + Au system at 
all energies considered here (150 AMeV to 2 AGeV) depending on the treat- 
ment of the binary collisions. The inclusion of a reaction plane conserving 
mechanism for the individual elastic particle-particle collisions has some ef- 
fect either on very small systems (cf. Fig. |5]) or when selecting particles with 
a low collision number in reactions between heavier systems. 



3.3 Mass distribution of resonances 

Starting again from the observation that the participating matter is the ori- 
gin of the transverse flow, another effect becomes worth pointing out: The 
amount of flow at high energies depends on the mass distribution of the res- 
onances, i.e. on the number of nucleon resonances taken into account in a 
transport model and being excited in the high density phase. 

Recent data on proton flow [38] indicate a decrease of sidewards flow 



above 1 AGeV incident energy following the well known logarithmic increase 
at low energies. Using standard potential parameterizations, both nonrela- 
tivistic 0, H, H] and relativistic [^|, |39j , this behavior cannot be understood 



within conventional transport models. In the latter the optical potential 
stays constant or even increases at high momenta and therefore the repul- 
sion generated from the momentum-dependent forces in a HIC gives rise to 
a significant contribution to the flow signal. However, since the nucleon- 
nucleus optical potential is only known up to 1 GeV experimentally ||40|| , it 



was recently proposed that this decrease of flow above 1 AGeV might indi- 
cate a decrease of the optical potential at high relative momenta or at high 
baryon density [[39] . 



In [11, 41] the transverse momentum of the baryons has been disentangled 
into a collisional part, a mean-field part and a part originating from the 
Fermi-motion of the particles, i.e. 

pt = p fl + p MF + p Fermi _ 

Depending on the EoS and rapidity interval considered, the relative contri- 
butions of these vary between 27, 73 and % for y ~ 0.5y pro and 37, and 
63 % for y > y pro for intermediate impact parameters. For more central colli- 
sions the Fermi-part contributes at minimum about 25 % over all rapidities. 
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'max 



However, since we deal with intermediate impact parameters b ~ 0.5 b r , 
only, which correspond to the multiplicity bins M3 and M4 where maximum 
transverse flow is produced, we neglect in the following the contribution of 
the Fermi-motion to the flow. 

Thus disentangling the flow signal into a collisional and a potential part, 
it turns out that ~50 % of the flow stems from the particle-particle collisions 
while roughly another 50 % are generated by the potential repulsion. Fig- 
ure P shows these two contributions for flow for the system Ni + Ni at b=4 



fm in comparison to the experimental data as compiled by Ref. |5S| using 



different EoS denoted by 'hard', 'medium' and 'soft'. Both the collisional 
"background" and the potential part rise up to 1 AGeV incident energy and 
remain constant above, whereas the data indicate a decrease above 1 AGeV. 
As seen from Figure |^, the Ni + Ni data are described reasonably by both 
a 'soft' and 'medium' EoS. 

Fig. [7], furthermore, shows in the upper part the flow in the cascade mode 
of the CBUU model, i.e. the collisional part, for the Ni + Ni system with 
all resonances up to M=1.95 GeV included (solid line) and without the res- 
onances above the A(1232) (dashed line). In the latter case all inelastic 
scattering strength is put into the NN — > iVA-channel, thus no stopping 
power from the inelastic collisions is lost. It is clearly seen that the exci- 
tation of higher nucleon resonances, though quite low in number, leads to 
an enhanced flow above 1 AGeV. In the lower part of Fig. [?] the results of 
the RBUU model 0, [§, |2|, g3[ g§, where only the A(1232) is included, 
is shown by the solid line. The flow results practically coincide with those 
from the CBUU model when including only the A resonance. In addition, 
we have calculated the flow within the RBUU model by shifting the A mass 
artificially to 1500 MeV (dashed line) in order to simulate effects from a 
heavy resonance. Again the flow increases compared to a calculation with 
the A(1232) mass above 1 AGeV and thus demonstrates the sensitivity of 
this observable with respect to the inelastic channels taken into account. 

We have analysed the origin of flow as a function of time within the various 
models and found that this effect is caused mainly by two reasons: (a) the 
Pauli blocking of final scattering states in the initial phase of the reaction, 
which is reduced in case of high mass excitations, and (b) a harder meson 
spectrum emitted by the fireball in case of higher resonances since these 
mesons partly are absorbed by the spectator matter which achieves a larger 
momentum transfer. Thus high mass resonances have to be included in any 
transport study that attempts to draw conclusions on the EoS or momentum 
dependence of the optical potential in comparison to experimental flow data 
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above bombarding energies of about 1 AGeV. 



4 Radial flow 



Radial flow has been discovered ||5| when analyzing the flow pattern of very 
central events of HIC. In contrast to transverse flow up to about 70 % of 
the incident energy (stored in the hot compressed fireball) is released as 
ordered radial expansion of the nuclear matter. Thus the hope is to extract 
information especially on the compressibility of the EoS via the magnitude 
of the radial flow. 

4.1 Temperature and flow energy 

Experimentally the radial flow is characterized or fitted in terms of the 



Siemens-Rasmussen formula 46 



d 3 N ^ P/T fsinha . ^ m , m 1 . , 

mm ~ p ■ ^ E,T {— ■ w + T ) - T - < 18 > 

with 7 = (1 — /3 2 )~ 1//2 and a = j(3p/T, where (3 denotes the flow velocity 
and T characterizes some temperature. We follow the same strategy in our 
CBUU calculations and apply a least square fit to the CBUU nucleon spectra 
using Eq. fll8|) . 

The results of the CBUU calculations for central Au + Au collisions are 
shown in Fig. |8| as a function of the bombarding energy in comparison to 



the data from []47|, |48|, [49| . We find that the 'temperature' T is systemati- 
cally underpredicted in all schemes investigated (soft and hard EoS, with and 
without momentum dependent forces), and that the flow velocities are not 
correctly reproduced, being too low at low energy and crossing the experi- 
mental data around 800 AMeV. This might, on the one side, indicate a strong 
binding from the potential, which gives not enough repulsion at high densi- 
ties and overcompensates the collisional pressure from the fireball. However, 
on the other side, the nucleon spectra resulting from the CBUU calculation 
show a strong non-thermal component at low incident energies and are thus 
in contradiction to the physical picture behind Eq. ( |18|) which assumes an 
isentropic expansion of a thermal equilibrated source. In a recent publication 



Ui~6 | we have investigated the degree of equilibration in a HIC as a function 
of the incident energy and the system mass and have found that even the 
most massive systems like Au + Au do not equilibrate at low energies. This 



is also in line with earlier findings [10 
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Thus we also evaluate the radial flow energy according to a prescription 
used in [51|] , 




which then gives the thermal energy as the difference to the total kinetic 
energy of the system 

1 A 



Etherm = ^ J2 [ V m i + Pi ~ m *J ~~ E fl™> " ( 20 ) 

This direct evaluation of Ef[ ow and Etherm or /5 and T, respectively, has the 
additional advantage of being statistically more stable than the fits to the 
nucleon spectra using Eq. (|18|). 

However, although Eq. fll^D obviously yields the radial motion of the 
system, the resulting flow values P ra diai differ strongly from the ones using 
Eq. as shown in Fig. |9] for central Au + Au collisions. Additionally, the 
results depend strongly on the evaluation time, since particle momentum and 
coordinate position become more and more aligned throughout the expansion 
phase. This is shown in Fig. |10] for a central Au + Au collision at 1 AGeV for 
the parameter T as well as the parameter (3fi ow . The changes of the results 
with time when using Eq. flUf) are much less pronounced. 

Consequently, one might ask for an "optimal" time for evaluating these 
quantities, at least when using Eqs. ( p!9| , p0|) . Looking at the time evolution 
of the central density, collision number and tt, A - abundancy throughout a 



central Au + Au collision at 1 AGeV (Fig. [11), we find that for this reaction 
the collision rate drops practically to zero at t ~ 25 fm/c. This is also the 
time when the average expansion velocities in transverse and longitudinal di- 



rection become approximately the same as shown in Fig. [L^ and the velocity 
profile reaches its final shape (Fig. |I^). Thus, inspite of further pion pro- 
duction via resonance decays, the reaction dynamics determining the flow is 
basically over at t ~ 25 fm/c and this time might be considered as a kind of 
hadronic freeze-out time. Nevertheless, the motion of the baryons becomes 
more ordered during the further expansion process as can be seen in Fig. [3~2| 
Consequently, when using Eq. (|19D Ef iow rises with time while E therm drops 
with time. However, already at the 'freeze-out time' t < 25 fm/c the results 
using Eqs. (|19l , |2~0"D and Eq. (|1|) differ quite substantially. 



Thus, in addition to Eq. (|18|) and Eq. ( PQ) for evaluating the 'temperature' 
T we have determined a 'temperature' via the N/ A - ratio. Since in |Tj| we 



have shown that at 1 AGeV at least isospin-equilibration is reached in a 
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Au + Au collision, we assume that also local thermal equilibrium is achieved 
and the N/ A - ratio might give a hint to the 'real' temperature of the system. 
The temperature T can be determined from 

Na 16 J dMf(M) M 3 / 2 e - E tvM)IT 

iVV ~ T m 3/2 e -B(p)/T ' ( ) 

where m,M denote the nucleon and A mass, respectively. The factor 16/4 
comes from spin-isosopin and 

2 M 2 r(M) 
^ ' ~ n (M 2 - M 2 R ) 2 + M 2 r 2 (M) (22) 

is the normalized A-mass distribution using 

w^-^feHS^)' ■ (23) 

as the momentum dependent A- width with £ 2 = (Mr — tjin — iti-k) 2 + rfj/4, 
Tr=150 MeV and p denoting the pion momentum in the Delta rest frame. 



Fig. |LJ shows the corresponding temperature in terms of the dotted line in 
comparison to the temperatures obtained using Eqs. (|T8| ) (full line) and (|20| ) 
(dashed line) as a function of time again for Au + Au at 1 AGeV. Though 
the N/ A-temperature shows a similar time dependence as the temperature 
obtained from Eq. (|20|), its actual values are closer to the temperatures ob- 
tained from the nucleon spectra over the whole time range. Especially at 
t=25 fm/c, where the collision rate drops to zero and the central density in 
the reaction volume gets below po, it is practically identical to the tempera- 
ture from Eq. fli"8D . Thus one might define t w 25 fm/c as the freeze-out time 
for this particular reaction. The nucleon spectra parameter T then conserves 
the N/A - 'temperature' at this time. 



5 Mass dependence 



As reported in Refs. 0, |57|, [58[ flow signals, both radial and transversal, are 
more pronounced for fragments with charge Z=2,3,- ■ ■ than for free protons 
or neutrons. This is usually attributed to i) the superposition of the ordered 
collective flow with the unordered thermal motion of the particles and ii) 
to the subsequent decay of excited heavier fragments during the expansion 
phase of the HIC into fragment plus nucleons: F'^ A+n ^ — > F(a) + n X N. 
Thus a direct comparison of the results from a single-particle model like 
CBUU, which describes the time evolution of the single-particle phase-space 
density, to the experimental data averaged over different fragments E7| seems 
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questionable. We thus study how far one can come in describing fragment 
flow data within the single-particle approach and then estimate the role of 
many-body correlations by the remaining deficiences. For transverse flow this 
approach has been used earlier [52, 53 in comparison to the experimental 
data from Refs. [||, |55|. 

We have applied a simple phase-space coalescence model to the final state 
distribution of the testparticles to form fragments out of the single-particle 
distribution. In this respect particles are decided to belong to a fragment, 
if their relative distance is less than r in coordinate space and smaller than 
Po in momentum space. While po was fixed to 250 MeV, which roughly 
corresponds to the average Fermi momentum of light nuclei, r was chosen 
to reproduce the experimentally measured number of free nucleons in the 
collision. 



The results for a central 1 AGeV Au+Au collision is shown in Fig. [15] for 
both the temperature T and flow velocity /3 ra diai derived from the respective 
spectra using Eq. (|18 ) in comparison to the results of Ref. [47 1 for proton 
spectra and averaged fits. In addition, the CBUU results from fitting the pure 
single-particle spectra are also shown in terms of the full squares. As can be 
seen from the upper picture, when applying the coalescence model to the final 
single-particle distribution, the temperature T for the remaining nucleons 
starts at the single-particle temperature (full square) and then increases with 
fragment mass towards the averaged experimental value. This increase is 
stronger for a hard equation of state (denoted by 'h' or 'hmd', respectively) as 
for a soft EoS (V and 'smd', respectively) due to a stronger collectivity of the 
nuclear motion. Averaging over the CBUU results for the different masses, 
the resulting mean temperature is in good agreement with the averaged ex- 
perimental measurement while the single-particle temperature is close to the 
experimental value for protons. 

The situation is similar for the flow velocity which drops with increasing 
mass, again roughly in line with the experimental findings of Ref. j|7|] . Also 
here the average over the CBUU results in the mass range of A = 1-4 leads 
to a flow value closer to the averaged experimental value. The different 
EoS employed (hard/soft, with and without momentum-dependent forces) 
practically lead to the same result for the flow velocity /3 ra< n a i- 

Again the results are quantitatively different when using Eq. (|19|) for 
evaluating the flow velocity. Whereas from the analysis of the spectra no 
difference between EoS of different compressibilities can be inferred, the flow 
energy derived using Eq. (plj) clearly separates a hard and soft EoS. Figure [16 



shows the results of the CBUU calculation for the 1 AGeV Au+Au collision 
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for fragment masses of 1 to 4 in comparison to the data from [4~7|. This 
comparison favors a hard EoS (full squares). 



For transverse flow we have compared the flow angle |58 



e, = t»-'f^, (24) 

V du z J 

with u x and u z denoting the velocity in transverse and longitudinal direction, 
respectively, to the data from Ref . . Since the statistics of the transverse 



momentum distributions of the fragments obtained by the coalescence model 
is poor, we have evaluated Eq. ( ^4[) for a given mass without separating for 
charge additionally. To be comparable to Ref. [pBfl , we estimate the mass 
number of the different charges measured by multiplying the given charge 
number with the average mass number for that charge. The result is shown 
in Figure [T7| for a Au+Au collision at 400 AMeV in terms of the full squares 



whereas the data are given in terms of the open circles. For the higher masses 
(A > 4) the agreement between the data and the CBUU calculation is very 
good, however, in the low mass region the flow angle is overestimated in the 
CBUU plus coalescence model quite significantly. 

In summary, the application of a simple coalescence model can provide 
results about the mass dependence of the radial flow parameters T and flradiai 
that are both correct in the functional dependence on mass and in their ab- 
solute magnitude. Here the important aspect is the suppression of random 
thermal motion, which can be provided by an arbitrary clustering algorithm. 
However, for transverse flow the situation is different. Obviously the ex- 
perimental proton momentum distribution is heavily distorted by massive 
fragment (or spectator) decay that is not adequately described in the single- 
particle model. This contrasts to the findings in an earlier publication |K 



which compared to PlasticBall data [54, 55 . 



6 Summary 

In this paper we have explored the dependence of transverse and radial flow 
signals on various model inputs and evaluation prescriptions using the CBUU 
model. 

For the transverse flow its origin could be traced back to the expand- 
ing participating matter. Starting from this observation we investigated the 
influence of different EoS, the conservation of angular momentum in the indi- 
vidual particle-particle collisions and of the mass distribution of the nucleonic 
resonances. 
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The most important point is that in general angular momentum conser- 
vation in the individual particle-particle collisions can be neglected for the 
overall description of a HIC. In contrast to that the mass distribution of the 
resonances included in the model - though the high mass resonances are quite 
low in number throughout a HIC - plays an important role for the descrip- 
tion of transverse flow above 1 AGeV. Since they change the collisional part 
of the signal substantially, they also influence the conclusions that one might 
draw from the experimental flow measurements on the nuclear potential and 
the EoS. 

For the radial flow we have concentrated on the difference between the 
results for the flow temperature T and flow velocity (3 when using different 
EoS and different evaluation prescriptions. The latter we see as an additional 
problem to the dependence of calculational results on nonphysical model 
inputs as recently pointed out also by Hartnack et al. ||. We showed that 
some kind of 'hadronic freeze out time' can be defined when looking at the 
collision rate, the isotropy of expansion and the flow velocity profiles and that 
at least at energies around 1 AGeV and heavy systems the thermodynamic 
temperature given by the A/A-ratio at this time is best reflected by the 
'temperature' of the nucleon spectra. 

In a third part we have investigated the fragment-mass dependence of the 
flow signals and have shown that even a simple coalescence model can provide 
reasonable agreement with the experimental findings. From this one might 
conclude that at least in central reactions of heavy systems (— > radial flow) 
we have a thermal scenario around 1 AGeV, while the difference between 
model and experiment for semicentral events (— > transverse flow) hints to- 
wards strong nonthermal contributions to the signals measured. 
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Figure captions 



Fig. [I]: Total energy per nucleoli (nucleon restmass subtracted) for a central 
Au + Au collision at 1 AGeV. Upper part: cascade mode; middle part: Mean- 
field propagation; lower part: full calculation including mean-field propaga- 
tion and collisions. 



Fig. ^ Coordinate space picture of a Au + Au collision at 1 AGeV at b=6 fm 
for different times. Initial stage (upper left), intermediate (upper right) and 
final stages (lower plots). The contour lines represent densities of 0.1, 0.5, 1, 
1.5 and 2 times po! the arrows indicate the direction and velocity of motion. 
In the lower plots the participating matter escaping from the reaction zone 
circumfloats the spectators. 

Fig. ^ Angular momentum per nucleon in a b=4 fm Ni + Ni collision at 150 
AMeV (upper part) and 2 AGeV (lower part). Shown are the cascade mode 
(dotted line), Vlasov mode (dashed line) and the full-BUU results (solid line). 

Fig. |3]: Angular momentum per nucleon in a b=4 fm Ni + Ni collision in the 
cascade mode. Upper part: Normal mode (dotted line) and with reaction 
plane conservation in the individual baryon-baryon collisions (solid line) at 
150 AMeV. Lower part: Normal mode (dotted line) and with relocation of 
the absorbing baryon to the CM in a meson-nucleon reaction (solid line) at 
2 AGeV. 



Fig. |5|: Azimuthal distribution of particles emitted into the forward hemi- 
sphere for systems of different size colliding with b ~ 0.6 x b max where b max 
is 2 x 1.125A 1 / 3 . The solid and dashed histograms are the CBUU results, the 
lines in the lower 3 panels are the appropriate fits in cos 2 (0) leading to the 
ratios of Table |2|. 

Fig. |(| Transverse flow for a Ni + Ni collision at b=4 fm as calculated within 
the CBUU model in cascade mode (dashed lower line) and for equations of 
state with different compressibilities K in comparison to data from EOS, 



Plastic Ball and FOPI as compiled by Ref. 38 



Fig. |7|: Transverse flow in the cascade mode from different transport models, 
CBUU (upper part) and RBUU ||4j |28| (lower part). The computations 
have been performed with different numbers of baryon resonances and reso- 



nance properties (see text). 



Fig. |8|: The radial flow velocity (lower part) and temperature (upper part) 
for central Au + Au collisions evaluated via Eq. (|18f) from the CBUU calcu- 
lations in comparison to the experimental data from Refs. |^7|, 48, The 



symbol V denotes a soft EoS without momentum dependent forces, 'h' a 
hard EoS and 'smd', 'mmd' and 'hmd' correspond to a soft, medium and 
hard momentum dependent EoS, respectively 

Fig. ^ The radial flow velocity for central Au + Au collisions evaluated ac- 
cording to Eq. ( |T8"D and Eq. ( |20| ) from the CBUU calculations as a function 
of the bombarding energy per nucleon E/A. 



Fig. 10: The radial flow velocity (upper part) and temperature (lower part) 
for a Au + Au collision at 1 AGeV evaluated according to Eq. (|j|) and 
Eqs. <MM) 

clS 8b function of the reaction time. 



Fig. [ll]: Time evolution of the central density, collision rate and pion num- 
ber throughout a central Au + Au collision at 1 AGeV. The thin vertical line 
marks t=25 fm/c. 



Fig. |12 : Time evolution of the total and radial mean velocity in transverse 



and longitudinal direction throughout a central Au + Au collision at 1 AGeV. 
The total velocity is given by the radial velocity by J2i 



Fig. |13|: The radial transverse velocity profile at different times of a central 
Au + Au collision at 1 AGeV. The final shape of the distribution is reached 
between t=22.5 and 25 fm/c. 

Fig. |14| : Temperatures obtained from the analysis of nucleon spectra using 
Eq. (|l^) (solid squares); direct evaluation according to Eqs. (|19|J20|) (open 
squares) and from the N/A ratio (dotted line). 



Fig. [15| : Fragment mass dependence of the parameters temperature T and 
flow velocity (3 from particle spectra using the CBUU plus coalescence model. 
Also shown are the results for the fit to the CBUU single-particle spectra (s.p., 
left) and the fits to the experimental spectra from Ref. [^7|| (exp., right). 



Fig. |16|: Flow energy for fragments of different mass evaluated using Eq. ([19]) 



for a central 1 AGeV Au + Au collision. The experimental data are taken 
from Ref. H7I . 



Fig. 17: Transverse flow angles for different mass fragments in a 400 AMeV 
Au + Au collision using the CBUU plus coalescence model in comparison to 



the data from Ref. [pq] . The horizontal lines indicate the flow angle satura- 
tion value. 
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